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O ■ ABSTRACT 

25 I We present the detection of nuclear stellar discs in the low-luminosity elliptical galaxies NGC 



o 
o 



4458 and NGC 4478, which are known to host a kinematically-decoupled core. Using archival 
HST imaging, and available absorption line-strength index data based on ground-based spec- 
troscopy, we investigate the photometric parameters and the properties of the stellar popula- 
tions of these central structures. Their scale length, h, and face-on central surface brightness, 
Qh' Hq, fit on fiQ-h relation for galaxy discs. For NGC 4458 these parameters are typical for nu- 

Q . clear discs, while the same quantities for NGC 4478 lie between those of nuclear discs and the 

• discs of discy ellipticals. We present Lick/IDS absorption line-strength measurements of H/3, 

^ I Mg 6, (Fc) along the major and minor axes of the galaxies. We model these data with simple 

i stellar populations that account for the a /Fc overabundance. The counter-rotating central disc 

of NGC 4458 is found to have similar properties to the decoupled cores of bright ellipticals. 
This galaxy has been found to be uniformly old despite being counter-rotating. In contrast, 
the cold central disc of NGC 4478 is younger, richer in metals and less overabundant than 
5— i ' the main body of the galaxy. This points to a prolonged star formation history, typical of an 

undisturbed disc-like, gas-rich (possibly pre-enriched) structure. 

Key words: galaxies: elliptical and lenticular, cD — galaxies: photometry — galaxies: kine- 
matics and dynamics — galaxies: abundances — galaxies: formation — galaxies: evolution 
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1 INTRODUCTION 

In recent years, the sub-arcsec resolution of the Hubble Space 
Telescope has allowed the study of galactic nuclei, unveiling the 
presence of distinct components such as small-scale stellar discs 
(see Pizzella et al. 2002 and references therein) and nuclear clus- 
ters (CaroUo et al. 1997; Boker et al. 2002; Graham & Guzman 
2003). To date surface-brightness distributions of nuclear stellar 
discs (hereinafter NSD) have been measured for only a few of SOs 
(Scorza &. van den Bosch 1998; van den Bosch et al. 1998; Ko- 
rmendy et al. 1996) and early-type spiral galaxies (Pizzella et al. 
2002). They have smaller scale lengths (~ 10-30 pc) and higher 
central surface brightnesses (~15-19 mag arcsec^^ in the V band) 
with respect to those of embedded stellar discs, which have been 
found in several elliptical and lenticular galaxies {e.g., Scorza & 
Bender 1995, hereinafter SB95). Although only few NSDs have 
been studied in detail, they may be quite a common structure in 



the central regions of spheroids. Indeed, nuclear discs of gas, dust 
and stars have been clearly detected in a large number of early-type 
galaxies {e.g., laffe et al. 1994; Lauer et al. 1995; van Dokkum 
& Franx 1995; Faber et al. 1997; Carollo et al. 1997; Tomita et 
al. 2000; Kormendy et al. 2001; Tran et al. 2001; Trujillo et al. 
2004), but the photometric parameters of such discs have not yet 
been derived. Furthermore, HST imaging in near-infrared band- 
passes provided indirect signatures of the presence of NSDs as 
discy isophotes in the nuclei of early-type galaxies (Ravindranath 
et al. 2001) or photometrically-distinct exponential components in 
bulges (Balcells et al. 2003, 2004). The existence of NSDs sug- 
gests that the continuity of the disc properties, with a smooth vari- 
ation of scale parameters from spirals to discy ellipticals along a 
sequence of decreasing disc-to-bulge ratio (Kormendy &. Bender 
1996), could be extended to nuclear scales (van den Bosch 1998). 
In this framework NSDs could provide important clues to the as- 
sembly scenario of their host galaxies. 
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In the current picture, the NSDs observed in bulges of SOs and 
spiral galaxies are believed to have formed either from a secular 
evolution of a nuclear bar [e.g., NGC 4570, van den Bosch et al. 
1998; Scorza & van den Bosch 1998; van den Bosch & Emsellem 
1998) or as the end result of a merging event {e.g., NGC 4698, 
Bertola et al. 1999; Pizzella et al. 2002). Each of these scenarios is 
likely to be correct for some but not for all the objects. In both 
of them the gas is efficiently directed toward the galaxy centre, 
where it first dissipates and settles onto an equilibrium plane and 
then forms into stars. Although the decoupled kinematics strongly 
suggest a later infall (e.g., NGC 4486A, Kormendy et al. 2001) 
or merger (e.g., Holley-Bockelmann Sl Richstone 2000), the fre- 
quently indistinguishable star formation history remains an enigma. 
Yet, these processes are expected to lead to distinct features in stel- 
lar populations of NSDs and surrounding spheroids. The measure- 
ment of the kinematics of some early-type galaxies (elliptical and 
SO galaxies) revealed evidence for kinematically decoupled or pe- 
culiar galaxy cores ( see Bertola & Corsini 1999 for a review), 
strongly indicative of formation by galaxy-galaxy merging (e.g., 
Kormendy 1984; see Mehlert et al. 1998) and references therein). 
Such kinematically-distinct cores have been associated with a cen- 
tral disc component (e.g., Rix & White 1992; Mehlert et al. 1998). 
However, triaxial galaxies with radius-dependent ellipsoids, can, in 
projection, mimic kinematically decoupled cores (see Arnold, de 
Zeeuw & Hunter 1994, and references therein). Spectroscopic ab- 
sorption line-strength indices (Faber et al. 1985; Trager et al. 1998) 
constrain the luminosity-weighted stellar population age, metal- 
licity and alpha-element overabundance by comparison with evo- 
lutionary population synthesis models (Worthey 1994; Maraston 
1998; Tantalo et al 1998; Vazdekis 1999; Thomas, Maraston & 
Bender 2003, hereafter TMB; Bruzual & Chariot 2003). Absorp- 
tion line-strength indices have been measured out to one effec- 
tive radii for elliptical galaxies (e.g., Davies et al. 1993; Gonzalez 
1993; Halliday 1999; Mehlert et al. 2000). Such measurements can 
distinguish between merging and dissipational formation scenar- 
ios by the constraint of metallicity gradients, dating of the latest 
episode of star formation and the measurement of the timescale 
over which the bulk of the star formation has taken place. Davies 
et al. (1993), Mehlert et al. (2003) and the analysis of literature 
data by Kobayashi & Arimoto (1999) have found elliptical galaxies 
to have gradients in metallicity shallower than predicted by dissi- 
pational collapse models. Mehlert et al. (1998) presented spectro- 
scopic line-strength data and HST archive imaging for two Coma 
cluster galaxies with kinematically-decoupled cores and found lit- 
tle evidence that the cores have experienced different star forma- 
tion histories from their parent galaxies. Using wide-field spec- 
troscopic maps of NGC 4365 acquired using the SAURON in- 
strument, Davies et al. (2001) presented evidence for a metal- 
rich core but constant age and a/Fe-enhancement for both the 
kinematically-distinct core and parent galaxy. Both studies are con- 
sistent with the core and the main galaxy having experienced simi- 
lar star formation histories. 

In this paper, we use archival HST imaging and ground-based 
spectroscopy to investigate the photometric parameters and stellar 
population diagnostics of the NSDs hosted by the two Virgo galax- 
ies NGC 4458 and NGC 4478. NGC 4458 and NGC 4478 are clas- 
sified as EO-1 (de Vaucouleurs et al. 1991, hereinafter RC3) and 
E2 (Sandage & Tammann 1981; RC3), respectively. They are both 
low-luminosity ellipticals with a M% ~ — 18 (RC3) at a distance 
of 12.6 Mpc (TuUy 1988, Ho = 100 km s"^ Mpc~^). They have a 
power-law central luminosity profile ( Ferrarese et al. 1994; Faber 
et al. 1997; Rest et al. 2001) and both galaxies show deviations 



from the r^''^ law in their outskirts (Michard 1985; Prugniel et al. 
1987; Caon et al. 1990; Peletier et al. 1990), which have been ex- 
plained as due to tidal interaction with NGC 4461 and NGC 4486. 
Recently, the stellar kinematics along the major and minor axes 
of NGC 4458 and NGC 4478 have been measured by Halliday et 
al. (2001), who found the signature of a kinematically-decoupled 
core in their inner ^ 5 arcsec. In particular, NGC 4458 has a clear 
counter-rotating core along the major axis, while NGC 4478 has a 
cold component detected along both the major and the minor axis. 

In Section |2| we perform the photometric analysis of the 
archival HST images of the nuclei of NGC 4458 and NGC 4478 in 
order to derive the surface brightness profiles and scale parameters 
of their NSDs. In Section|3|we present the absorption line-strength 
indices as a function of radius of the major and minor axes of both 
galaxies from the thesis research of Halliday (1999). In Section|4] 
we combine photometric and spectroscopic results to assess age, 
metallicity and overabundances of stellar populations of the two 
NSDs and surrounding spheroids. This allowed us to suggest pos- 
sible formation and evolution scenarios for NGC 4458 and NGC 
4478. We summarize our conclusions in Section|5| 



2 PHOTOMETRIC PARAMETERS 
2.1 Data reduction 

We retrieved Wide Field Planetary Camera 2 (WFPC2) images of 
NGC 4458 and NGC 4478 from the HST archive. Data for the filter 
F814W were selected as a compromise between obtaining data in 
identical bandpasses for both galaxies and minimizing the effects 
of dust on photometric measurements. Total exposure times were 
1 120 s for NGC 4458 (5 exposures. Prog. Id. 55 12, PI. S. M. Faber) 
and 1600 s for NGC 4478 (2 exposures. Prog. Id. 6587, PI. D. Rich- 
stone), respectively. All exposures were taken with the telescope 
guiding in fine lock, which typically gave an rms tracking error 
of 0.003 arcsec. We focused our attention on the Planetary Camera 
chip (PC) where the nucleus of both galaxies was centred. This con- 
sists of 800 X 800 pixels of 0.0455 x 0.0455 arcsec^ each, yielding 
a field of view of about 36 x 36 arcsec^ . The images were calibrated 
using the standard reduction pipeline maintained by the Space Tele- 
scope Science Institute. Reduction steps include bias subtraction, 
dark current subtraction, and flat-fielding and are described in detail 
in Holtzman et al. (1995a). Subsequent reduction was completed 
using standard tasks in the STSDAS package of IRAF^. Bad pixels 
were corrected by means of a linear one-dimensional interpolation 
using the data quality files and the WFIXUP task. Dilferent images 
of the same target were aligned and combined using IMSHIFT and 
knowledge of the offset shifts. Cosmic ray events were removed 
using the task CRRE J. The cosmic-ray removal and bad pixel cor- 
rection were checked by inspection of the residual images between 
the cleaned and combined image and each of the original frames. 
Residual cosmic rays and bad pixels in the PC were corrected by 
manually editing the combined image with IMEDIT. The sky level 
(~ 1 count pixel" ^) was determined from regions free of sources 
in the Wide Field chips and subtracted from the PC frame after ap- 
propriate scaling. 



^ IRAF is distributed by NOAO, which is operated by AURA Inc., under 
contract with the National Science Foundation 
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Figure 1. Left panels: WFPC2/F8 14W images of NGC 4458 and NGC 4478 
(for which nuclear discs were detected), and NGC 4539 (for which a nuclear 
disc was not detected). The size of the plotted region is 19.3 X 19.3 arcsec^ . 
The orientation is specified by the arrow indicating north and the segment 
indicating east in the lower right comer of each panel. Middle and right 
panels: Unsharp masking of the WFPC2/F814W images obtained with a = 
2 and 6 pixels, respectively. Sizes and orientations are as in the left-hand 
panels. 



2.2 Detection of nuclear discs 

To test for the presence of a NSD in NGC 4458 and NGC 4478, we 
constructed the unsharp-masked image of the PC frame using an 
identical procedure to Pizzella et al. (2002). We divided each image 
by itself after convolution by a circular Gaussian of width a = 2 
and 6 pixels, corresponding to 0.09 arcsec and 0.27 arcsec, respec- 
tively (Figure0. This procedure enhanced any surface-brightness 
fluctuation and non-circular structure extending over a spatial re- 
gion comparable to the a of the smoothing Gaussian. Two values 
of a were adopted to attempt the identification of structures with 
different scale lengths. 

A highly elongated structure is clearly visible in the nucleus 
of each galaxies. These structures are not artifacts of the unsharp- 
masking procedure. In FigureQwe show as a counter-example the 
unsharp-masked image of NGC 4539 which has been obtained us- 
ing an identical method to that for NGC 4458 and NGC 4478; we 
clearly do not detect an elongated disc-like structure. Moreover the 
nuclear structures of NGC 4458 and NGC 4478 are associated with 
a central increase in ellipticity measured by performing an isopho- 
tal analysis using the I RAF task ELLIPSE (Figure|2j. 

We analyzed the isophotal profiles of both galaxies by first 
masking foreground stars and then fitting ellipses to the isophotes. 
We allowed the centres of the ellipses to vary, to test whether the 
galaxies were disturbed. Within the errors of the fits, we found no 
evidence of variations in the fitted centre. The ellipse fitting was 
repeated with the ellipse centres fixed. The resulting azimuthally 
averaged surface brightness, ellipticity, position angle, the fourth 
(04), and sixth (ae) cosine Fourier coefficients profiles are pre- 
sented in Figure|2| In the innermost ~ 2 arcsec we measured pos- 
itive values of the 04, and ag Fourier coefficients, which describe 
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Figure 2. (a) Radial profiles of surface brightness of NGC 4458 {left panel) 
and NGC 4478 (right panel) after deconvolution (open squares) are com- 
pared to those of their nuclear discs (thick lines), (b) Radial profiles of the 
fraction of total luminosity contributed by the nuclear disc of NGC 4458 
(left panel) and NGC 4478 (right panel). Radial profiles of (cj ellipticity, (d) 
position angle, fej 04 and (f) ag Fourier coefficients before (open squares) 
and after (filled circles) the subtraction of the best-fitting model for the nu- 
clear disc of NGC 4458 (left panels) and NGC 4478 (right panels). 



the discy deviation of the isophotes from pure ellipses (Jedrzejew- 
ski 1987). These photometric features have been also observed in 
the HST/F555W image for NGC 4458 (Lauer et al. 1995; Trujillo 
et al. 2004) and in the HST/F702W image of NGC 4478 (van den 
Bosch et al. 1994; Rest et al. 2001 ; Trujillo et al. 2004) and confirm 
the presence of a NSD in both galaxies. 



2.3 Photometric decomposition 

2.3.1 Photometric parameters of the miclear discs with the 
Scorza & Bender method 

After establishing the existence of NSDs, we measured their pho- 
tometric properties using the method described by SB95. This 
method is based on the assumption that isophotal disciness is the 
result of the superimposition of a spheroidal component (which is 
either a host elliptical galaxy or a bulge component) and an in- 
clined disc. The two components are assumed to have both per- 
fectly elliptical isophotes with constant but different ellipticities. 
When adopting this technique to study the innermost regions of 
galaxies it is essential to restore the images from the effects of 
the HST point spread function (PSF) in order to properly derive 
the photometric parameters of the NSDs, as shown by Scorza & 
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van den Bosch (1998). Such deconvolution was performed through 
the Richardson-Lucy method by means of the IRAF task LUCY. 
Although susceptible to noise amplification, this algorithm was 
shown by van den Bosch et al. (1998) to lead to a restored surface- 
brightness distribution comparable to that obtained using a multi- 
Gaussian representation (Monnet et al. 1992; Cappellari 2002). We 
believe that the results obtained by van den Bosch et al. (1998) 
are directly applicable to our case: we are dealing with images ob- 
tained with similar or longer integration times, galaxies with less 
steep surface-brightness profiles and NSDs with equal or larger 
scale lengths. We decided to deconvolve the images with a num- 
ber of iterations between 3 and 6. A larger number of iterations 
does not affect the result of the decomposition but does amplify 
the noise. For each given image and nucleus position on the PC 
we adopted a model PSF calculated using the TINY-TIM package 
(Krist & Hook 1999). No correction for telescope jitter was nec- 
essary. The SB95 method consists of the iterative subtraction of a 
thin disc model characterized by an exponential surface-brightness 
profile with central surface brightness Jo and radial scale length h, 
and by an axial ratio b/a. The disc parameters are adjusted until the 
departures from perfect ellipses are minimized (i.e., 04 and are 
close to zero). For each disc model we obtained the disc-free image 
of the galaxy by subtracting the disc model from the galaxy image. 
We performed an isophotal analysis on the disc-free image using 
the IRAF task ELLIPSE. We defined 

where a^^disc-free{i) ± o-{i) is the value of the 04 Fourier coef- 
ficient measured for the i— th isophote in the disc-free image, and 
A'' is the number of fitted isophotes. We assumed a{i) = 0.01 
as typical error on a^, disc- free for all the isophotes in the region 
of the NSD. We defined the reduced as xl = xV(^ - M), 
where Af = 3 is the number of free parameters, namely Jo, h, and 
b/a. According to our technique for photometric decomposition the 
minimum value of xt = x1,min corresponds to the best-fitting 
model of the NSD. We found' x^,™,„ = 0.43 for NGC 4458 and 
Xl,m,n = 0.34 for NGC 4478. We determined Ax2 = xl~xl,,mn 
and derived its confidence levels under the assumption that the er- 
rors are normally distributed (Press et al. 1992). The resulting con- 
tour plots of xt are shown in Figure|3|and they allowed us to derive 
the best-fitting values of Jo, h and b/a and their la errors. 

We derived the Johnson V— band central surface brightness 
of the NSDs from the HST VEGAMAG system (Holtzman et al. 
1995b) using the STSDAS task SYNPHOT. Since this correction 
depends on the galaxy spectral energy distribution, it has been cal- 
culated using the spectral template for elliptical galaxies by Kinney 
et al. (1996). The resulting relation is V — mF8i4w = 1.33. We 
derived the scale length of the NSDs assuming a distance of 12.6 
Mpc for both galaxies. The inclination of NSDs has been calculated 
as i = arccos {b/a). The resulting values of central surface bright- 
ness in the V band, scale length and inclination for both NSDs are 
listed in TableQ 

The disciness of NGC 4458 and NGC 4478 disappears after 
the subtraction of the best-fitting model of the NSD as it results 
from the elliptical shape of their isophotes as shown by the pho- 
tometric profiles plotted in Figure |2| We tested that the result of 
the photometric decomposition is independent of position on the 
PC chip where the PSF is built. The resulting fraction of total lu- 
minosity contributed by the nuclear disc of NGC 4458 and NGC 
4478 is plotted as a function of radius in Figure|2| 



Table 1. Photometric parameters of the nuclear di.scs 



Name ii^y h i L^„^ 

[mag arcsec~2] [pc] [°] [Lq.v] 

NGC 4458 14.74^°;^^ H-Ola.s 83-1^8.8 4.8 ■ 10«^ 

NGC 4478 le.eilg g* 40.5+5g ° ''7.9t",;g 1.9-10'' 



2.3.2 Photometric parameters of the nuclear discs with the 
Seifert & Scorza method 

In applying the decomposition method by SB95 we assumed that 
the surface-brightness profile of the NSDs was exponential. To ver- 
ify this assumption, we derived the surface-brightness profile of 
the NSDs using the alternative method of Seifert & Scorza (1996, 
hereinafter SS96). This decomposition method was developed to 
identify embedded discs in ellipticals and bulges without any a pri- 
ori assumptions of the parametric law for their surface-brightness 
profile. In the deconvolved image of the galaxy we masked the 
region of the NSD with a double-cone mask, and fit ellipses to 
the isophotes in the remaining part of image using the IRAF task 
ELLIPSE to derive the photometric properties of the host galaxy. 
A double-cone mask centred on the galaxy major axis is an optimal 
choice for the highly-inclined NSDs of NGC 4458 and NGC 4478 
since we expect that the isophotes close to the galaxy minor axis 
will be only slightly contaminated by light from the NSDs. The 
very central portion of the galaxy image was not masked to allow 
meaningful isophotes to be fitted. These masks correspond to the 
central 3x3 pixel (0.14 x 0.14 arcsec^) of NGC 4458 and 5 x 5 
pixel (0.23 x 0.23 arcsec^) of NGC 4478. The opening angle of the 
mask (~ 90°) was adjusted to minimize the influence of the NSD 
on the fitted parameters of the surrounding galaxy bulge compo- 
nent, but to ensure that a sufficient number of data points for the 
model derivation remained. A model was constructed to have per- 
fectly elliptical isophotes with the same surface brightness, ellip- 
ticity and position angle radial profiles which we obtained from the 
isophotal fit. We subtracted the galaxy model from the original im- 
age to obtain a residual image which shows an elongated structure, 
namely the NSD, responsible for the higher-order deviations from 
perfectly elliptical isophotes. We extracted the surface-brightness 
radial profile along the major axis of the residual disc which was 
close to being exponential. Scale lengths are consistent within their 
errors to those measured using the SB95 method, while the central 
surface brightnesses are somewhat fainter. This discrepancy in the 
surface brightnesses derived for using the SB95 and SS96 meth- 
ods is due to the actual shape of the double-cone mask where we 
excluded the central region in order to allow meaningful isophotes 
to be fitted. In this way the observed central surface brightness is 
assigned to the host galaxy and the central surface brightness of the 
NSD is therefore underestimated. 



2.4 Structural parameters of nuclear discs 

It is interesting to compare the scale lengths given in Table Q 
(h ^ 0.2" for NGC 4458 and h « 0.7" for NGC 4478) with the 
"break" radii found in the major and minor axis kinematics of the 
two galaxies (Halliday et al. 2001). Counter-rotation of about (30 
km s~^) is observed in the inner 4 arcsec of the major axis of NGC 
4458. In the case of NGC 4478 the velocity dispersion drops in the 
central 2 arcsec, where the rotational velocity along the major axis 
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Figure 3. Contour plots of xi for NGC 4458 (upper panels) and NGC 4478 (lower panels) as a function of the central surface brightness, Iq, scale length, h, 
and axial ratio, b/a, of the nuclear disc models. Solid contours define the formal 68.3%, 95.4%, and 99.7% confidence levels. Intermediate confidence levels 
are given with dashed contours. Filled circles indicate actual model calculations. Bi-cubic interpolation is used to calculate at intermediate points. Crosses 
correspond to the best-fitting photometric parameters of the nuclear disc. 



drops by « 20 km from 50 km at 2 arcsec to 30 km 
at 3 arcsec from the centre to increase again at larger radii. 

Clearly, the spatial resolution of the kinematic data (see Sec- 
tion |3| for details) is smearing the signal of the NSDs. Figure 4 
shows the position of the NSDs of NGC 4458 and NGC 4478 in 
the disc Ho-h diagram (adapted from Pizzella et al. 2002; but see 
also van den Bosch 1998, and Graham 2001). The face-on central 
surface brightness fig have been derived from the observed one no 
by applying the following inclination correction 

1^0 = MO —2.5 log(cosj) (2) 

and without taking into account any correction for extinction. The 
(small) scale length and (high) central surface brightness of the disc 
of NGC 4458 are typical for NSDs. In contrast, the (relatively large) 
scale length and (relatively low) central surface brightness of the 
disc of NGC 4478 are in between the typical values found for discy 
ellipticals and NSDs. No correlations are found with global prop- 
erties (e.g., luminosity, scale-length, disc-to-bulge ratio) of the host 
galaxies. 

Finally, on a speculative line of thought without presumption 
of rigourness, it is interesting to correlate the order of magnitude 
rotational velocities estimated above (30 km s^^ for NGC 4458 
and 50 km for NGC 4478) with the estimated magnitudes of 
the discs. We find that discs following the TuUy-Fisher relation (as 
in Haynes et al. 1999) and having the same rotational velocities, 
are expected to be ~ 3 mag brighter than the NSDs of NGC 4458 
and NGC 4478. In other words, the NSDs of NGC 4458 and NGC 
4478 rotate faster than discs of the same luminosity that follow the 
Tully-Fisher relation. Taken at face value, this might indicate that 
they behave like test particles embedded in the dominant gravita- 



tional potential generated by the main (bulge) body of the galaxies. 
However, the situation is complicated by the low spatial resolution 
of the kinematic data. Using the total masses of the discs derived in 
Section|4| we estimate that maximum rotational velocities of the or- 
der of ~ 30 — 50 km are generated by the NSDs alone around 
0.6-2 arcsec. Clearly, higher resolution kinematic data and proper 
modelling (including a decomposition of the line-of-sight velocity 
distributions and seeing convolution) are needed to pin down the 
issue. 



3 LINE-STRENGTH INDICES 

Absorption line-strength indices are taken from a study of the stel- 
lar populations and kinematics of low-luminosity elliptical galax- 
ies, in the Virgo cluster and nearby groups (Halliday 1999; Hal- 
liday et al. 2001). We provide a brief summary of our spectro- 
scopic observations. Long-slit spectroscopic data were obtained us- 
ing the Blue Channel Spectrograph at the Multiple Mirror Tele- 
scope, Arizona (USA). A Loral 3kxlk CCD chip with a pixel 
size of 15/imxl5/im, a grating of 1200 slits mm~^ and a slit 
width of 1 arcsec were used. The spectral resolution was ~1.5 A 
FWHM, wavelength range 4555 - 6045 A and spatial scale 0.6 arc- 
sec pixel"^. Seeing varied between 0.5 and 0.8 arcsec during ob- 
servations. The reader is referred to Halliday (1999) for details of 
the spectroscopic data reduction and Halliday et al. (2001) for a 
description of the kinematical analyses. 
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Figure 4. Disc ^l'^-h diagram adapted from Pizzella et al. (2002). 
Open circles refer to liigli surface-briglitness spirals, triangles to 
low surface-brightness spirals, stars to SO's, filled circles to discy 
ellipticals. Large squares correspond to NSDs in elliptical galaxies 
(open symbols) and in SO's and spiral galaxies (filled symbols), re- 
spectively. 



3.1 Measurement of absorption line-strength indices 

Line-strength indices were measured as a function of radius to 
approximately one photometric effective radius along the major 
and minor axes of each galaxy. One-dimensional spectra were ex- 
tracted for different luminosity-weighted radii of each axis and 
Une-strength indices were measured for each spectrum using an al- 
gorithm developed by G. Baggley and revised by H. Kuntschner 
(Halliday 1999). The Lick/IDS system line-strength index defini- 
tions of Trager et al. (1998) were adopted. A fuller description of 
the algorithm is given below. 

The principal input was a spectroscopic CCD frame for which 
a basic data reduction had been completed. Before extracting one- 
dimensional spectra the galaxy spectrum was checked to be aligned 
with the rows of the CCD frame. Using the I RAF tasks APFIND 
and APTRACE the position of the galaxy centre was traced along 
the dispersion axis; this step was performed to an accuracy of 
~0.03 pixels rms. The pixels were resampled such that the galaxy 
spectrum was centred on the same row pixel position as a function 
of dispersion. 

One-dimensional spectra were extracted by summing rows of 
the galaxy frame until a signal-to-noise ratio (hereinafter S/N) of 40 
per A (for the wavelength range ~5 100-5300 A) was attained. The 
dispersion axes of extracted spectra were rebinned to logarithmic 
intervals of wavelength. Redshifts were measured using the I RAF 
task FXCOR, by cross-correlating each spectrum with a range of 
one-dimensional de-redshifted stellar spectra. The mean redshift 
determined for all templates was used to de-redshift each galaxy 
spectrum. H/?, Mg b, Fe5270 and Fe5335 line-strength indices 
were measured for all spectra. Errors for these raw index measure- 
ments were calculated based on Poissonian noise, the typical sky 
level counts subtracted during sky subtraction, and the values of 
read-out noise and gain of the CCD. 

Our line-strength index measurements were transformed to the 
Lick/IDS system. This procedure involved two steps: (i) correc- 



tion to zero galaxy velocity dispersion and (ii) correction to the 
Lick/IDS spectral resolution (^8.6 A). 

To correct for the effects of galaxy velocity dispersion, a, 
different stellar spectra were smoothed by Gaussians correspond- 
ing to measurements of cr between km s~^ and 300 km s~^ 
in intervals of 20 km s^^. Smoothing was performed using the 
I RAF task GAU S S . Correction factors for different values of a were 
calculated for each index by comparing Une index measurements 
for each smoothed spectrum with measurements for the original 
unsmoothed stellar spectrum. For the atomic indices H/?, Mg b, 
Fe5270 and Fe5335 correction factors were defined to be the ra- 
tio of measurements for the original unbroadened stellar spectrum, 
lorig, and that for the spectrum broadened to a particular velocity 
dispersion, i.e., 

C{a) = (3) 

Mg b, Fe5270 and Fe5335 Une-strength index measurements were 
corrected for the effects of a using the calibrations for 12 stellar 
observations; for H,5 data 8 stellar observations wore used. Linear 
relations were interpolated between the mean correction factor for 
each a. The correction factor for an arbitrary measurement of cr was 
determined using the linear relation for the a range bracketing the 
a measurement. The measurement of cr as a function of radius was 
presented for both galaxies in Halliday et al. (2001). We adopt the 
cr measurements obtained for a S/N ~ 60 per A. Measurements 
of cr with S/N = 30 per A were adopted at large radii where 
uncertainties were larger. 

Absorption line-strength index measurements and measure- 
ments of cr from Halliday et al. (2001) were measured in general for 
different galactic radii. Linear interpolations were made between 
neighboring measurements of cr. A cr measurement for a radius at 
which line-strength indices were measured, was found by linear 
interpolation from bracketing data. Depending on the appropriate 
side of the galaxy, if the radius was either less than or exceeded the 
radius of the most reUable measurement of a, the Unear interpo- 
lation determined for the closest radial interval was used to deter- 
mine a. The errors of line-strength index measurements corrected 
for galaxy velocity dispersion were calculated to be, 

Slatormc(a) = ^/(Slorrg ■ C(cr))2 (^J^^^.^^ . ^^(cr))^ (4) 

where Slorig is the error in the raw measurement calculated as out- 
lined above and 6C{a) is the addition in quadrature of the standard 
deviations of correction factors for measurements of cr bracketing 
the cr measurement. 

Corrections for the difference between our spectral resolution 
and the resolution of the Lick/IDS system were applied. A simu- 
lated Lick aperture was extracted for the major axis spectrum of 
each galaxy studied by Halliday (1999) present in the "pristine" 
IDS sample of Trager (1997). Spectra were resampled such that 
the galaxy centre was aligned with the same pixel row position as 
a function of wavelength. For our spectroscopic slit width of 1 arc- 
sec and spatial scale of 0.6 arcsec pixel" ^ a one-dimensional spec- 
trum corresponding to a Lick/IDS aperture was created by sum- 
ming the central 9 rows (i.e., the central ^ 5.4 arcsec of each 
galaxy spectrum). Line-strength indices were measured using the 
Lick/IDS definitions of Trager et al. (1998). The galaxy veloc- 
ity dispersion was measured using the IRAF task FXCOR. Line- 
strength indices were corrected for the effects of a using the cali- 
brations discussed above. Direct comparisons were made with the 
fully-corrected Lick/IDS 1.4 x 4.0 arcsec^ aperture measurements 
of the "pristine" IDS sample of Trager (1997); we considered 
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the Lick/IDS H/3 measurements corrected for the effects of a but 
uncorrected for H/3 emission. Relation fitting was attempted be- 
tween the simulated Lick/IDS aperture measurements corrected for 
a and the Lick/IDS measurements of Trager ( 1 997), as a function of 
our a-corrected, simulated Lick/IDS aperture measurements. After 
testing various routines, the median value of offset was found to 
provide a good fit for all line-strength indices {i.e., for H/3, 0.272; 
Mg2, -0.048; Mgb, 0.051; Fe5270, 0.544; Fe5335, 0.697). The cr- 
corrected measurements were established on to the Lick/IDS sys- 
tem by subtracting these median values for all index measurements 
excluding Mg b for which no zero correction was required due to 
the large scatter of offset values around the median. 



3.2 Absorption line-strength gradients 

Absorption line-strength index measurements for H/3, Mg h, Mg2 
and (Fe) = (Fe5270 + Fe5335)/2 are shown as a function of ra- 
dius for the major and minor axes of NGC 4458 and NGC 4478 in 
Figure |5| Line-strength measurements are presented as a function 
of the normalized radius R/R*, where _R* = R^j \J\ — e for the 
major axis and R* — R^ ■ \/l — e for the minor axis, respectively. 
Re is the circularized effective radius defined as Re = Ae/2 with 
Ae the diameter of the circle enclosing one-half of the B— band 
total flux (RC3), and e is the mean ellipticity of the galaxy. We 
assume Re — 26.1 arcsec (RC3, corresponding to 1.6 kpc at the 
adopted distance) and e = 0.00 (Lauer et al. 1995) for NGC 4458, 
Re = 13.4 arcsec (RC3, corresponding to 0.8 kpc at the adopted 
distance) and e = 0.17 (Peletier et al. 1990) for NGC 4478. The 
different central values measured along major and minor axis for 
the same index give an estimate of the true measurement errors. 
Measurements for Mg2 are shown for completeness and the cali- 
bration of this molecular index to the Lick/IDS system is described 
in Halliday (1999). 

For NGC 4478, measurements of H/3, Mg fe, Mg2, Fe5270, 
Fe5335 and (Fe) were derived previously for the galaxy centre by 
Gonzalez (1993) and Trager et al. (1998), and as a function of ra- 
dius of the major axis by Peletier (1989) and Gorgas et al. (1990). 
In Figure |6| we compare all previous measurements with our mea- 
surements presented here. 

Good agreement within all margins of uncertainty is found be- 
tween our data and the measurements of Peletier (1989). Measure- 
ments and errors for Fe5270 were not tabulated by Peletier (1989) 
and have been derived using the corresponding measurements of 
Fe5335 and (Fe). Agreement within the errors is found between 
our measurements and those of Gorgas et al. (1990) for H/3 and Mg 
fe but not for Fe5270, Fe5335, and (Fe). A single bad pixel column 
was interpolated for a region corresponding to the red passband of 
Fe5270 of our spectrum of NGC 4478 but this is unlikely to be the 
cause of the discrepancy here since no evidence of unsuccessful 
bad pixel correction was found. Since Fe5270, Fe5335, and (Fe) 
are critical indices in the study of stellar populations, it has to be 
noted that our results are based on our data. At the galaxy centre 
there is a small difference between our value of the Mg2 index and 
that of Gorgas et al. (1990), while agreement within the errors is 
found for data at all other radii. The central measurements of the 
line-strength indices of Gonzalez (1993) and Trager et al. (1998) 
are in good agreement with our measurements. 
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Figure 6. Comparison of our line-strength index data for NGC 4478 mea- 
sured as a function of radius of the major axis (shown by filled triangles 
symbols), the data Gorgas et al. (1990, open squares) and Peletier (1989, 
asterisks). The central measurements of Gonzalez (1993, pentagon) and 
Trager et al. (1998, hexagon) are shown for representative radii. Error bars 
are provided for data at the innermost and outermost radii and the central 
measurements. 



4 AGES, METALLICITIES AND 
Q/Fe OVERABUNDANCES 

In the following we indicate the combined Magnesium-Iron index 
with 



[MgFe]' = ^Mg 6(0.72 x Fe5270 + 0.28 x Fe5335) 



(5) 



newly defined by TMB. This index is completely independent of 
the [a/Fe] overabundance and hence serves best as a metallicity 
tracer. 

We show in the top panels of Figure the line strengths of 
the H/3 and [MgFe]' indices at different positions in the galaxies 
NGC 4458 and NGC 4478 with the grid of models of TMB. No 
clear trend is observed for NGC 4458. NGC 4478 shows a negative 
gradient in both the indices H/3 and [MgFe]'. We show in bottom 
panels of Figure|7|the line strength of the Mg 6 and (Fe) indices at 
different positions in the galaxies NGC 4458 and NGC 4478. For 
both galaxies, both indices (Fe) and Mg b decrease outward on the 
mean. 

Using stellar population models with variable element abun- 
dance ratios from TMB we derive average ages, metallicities and 
[a/Fe] ratios as a function of the position in the galaxies, follow- 
ing the procedure of Mehlert et al. (2003). Figure [8| translates the 
trends observed in FigureQas follows. NGC 4458 is uniformly old 
{t ~ 15 Gyr), has overall low metal content ([Z/H] ~ —0.2) and 
overabundance ([a/Fe] ^ 0.3). For NGC 4478 instead the cen- 
tral regions {R < O.lRe) are younger (t ~ 6 Gyr), more metal 
rich ([Z/H] ^ 0.35) and less overabundant ([a/Fe] ^ 0.2) than 
the outer regions, where old ages (t ~ 15 Gyr), low metallicities 
([Z/H] ~ 0.1) and high overabundances ([a/Fe] ~ 0.3) are ob- 
served on both the minor and the major axis. The size of the error 
bars shows that the decreases in age and overabundance in the cen- 
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Figure 5. The radial profiles of the line-strength indices measured along the major and minor axis of NGC 4458 {left panels) and NGC 4478 {right panels). 
Squares refer to the value of line-strength indices measured at the galaxy centre. 



tral region are significant. Finally, note that within the errors most 
of the ages shown in Fig.|8|are compatible with the age of Universe 
favoured by WMAP (Spergel et al. 2003). In the following we will 
interpret the radial variations described above as the result of the 
superposition of the central discs detected photometrically in Sec- 
tion l2.2l and the main body of the galaxy. Moreover, we will assume 
that the main body of the galaxy does not have a [a /Fe] gradient, in 
agreement with the results of Mehlert et al. (2003). NGC 4458 has 
properties similar to giant ellipticals with decoupled cores. Proto- 
types of decoupled cores in giant ellipticals are as old as and metal 
richer than the rest of the galaxy, with high and approximately con- 
stant overabundance, as in the case of NGC 4816 and IC 4051 
(Mehlert et al. 1998), or NGC 4365 (Surma & Bender 1995; Davies 
et al. 2001). In NGC 4458 the inner disc and the main body of the 
galaxy appear to have the same stellar populations, with a metal- 
licity gradient for the main body of the galaxy slightly weaker than 
the cases discussed above. By contrast, the case of NGC 4478 is 
intriguing. The inner, cold disc is younger, richer in metals and less 
overabundant than the main body of the galaxy. We tested this con- 
clusion with a composite stellar population model, summing 10% 
of a young {t — 1.65 Gyr) and metal rich (Z = 3Zq) component 
with solar a/Fe ratios, and 90% of an old {t = 12 Gyr) metal rich 
{Z — 2Zq) overabundant ([a/Fe] = 0.4) population, and recover- 
ing the observed central values of the line indices. Given the errors 
of the available spectroscopic data and their spatial resolution, we 
refrained from a detailed study of the allowed range of alternative 
solutions, that only spectroscopy at HST resolution can pin down. 

Of course, higher than solar a/Fe ratios for the disk are al- 
lowed if we relax the assumption that the a/Fe ratio of the main 
body of NGC 4478 is constant. Using these age and metallicity es- 
timates and the models of Maraston (1998), we compute that the 
mass-to-light ratios of the inner discs are M / Lv = 7.2 M0/L0 
for NGC 4458 (with a range from 5 to 10 judged from the range of 
ages and metallicities shows in Fig.|8j and M / Lv = 1.3 M© /L© 
for NGC 4478. In this case, given the limitations of the spectro- 
scopic data discussed above, the uncertainties are probably as large 



as a factor 2. Combined with the luminosities given in Table □ 
this gives masses of M = 3.4 x 10^ M© and M = 2.5 x 10^ 
M0, respectively. The kinematic and photometric properties dis- 
cussed above argue for the presence of (cold) discs at the centres of 
NGC 4458 and NGC 4478, and therefore for a formation scenario 
through the dissipational collapse of a gas rich object. The observed 
kinematic signature of counter-rotation in NGC 4458 argues in ad- 
dition for an external origin with decoupled angular momentum. 
The absence of gradients in the stellar population of NGC 4458 in- 
dicates that the formation of the inner cold structure happened at 
the same time as the main body of the galaxy. The younger age 
and low overabundance of the central structure of NGC 4478 is in- 
dicative of a prolonged star formation history (Thomas, Greggio & 
Bender 1999), typical of an undisturbed disc-like, gas-rich (possi- 
bly pre-enriched) structure. This is consistent with the absence of 
counter-rotation. 



5 CONCLUSIONS 

We have investigated the properties of the central regions of the 
low-luminosity ellipticals NGC 4458 and NGC 4478, finding that: 

(i) HST archive images of the two galaxies show the presence of 
central discs. They have scale lengths and face-on central surface 
brightnesses which fit on ^0-/1 relation for galaxy discs. For NGC 
4458 these parameters are typical of NSDs, while for NGC 4478 
they are intermediate between the values for NSDs and the discs of 
discy ellipticals. 

(ii) Measurements presented by Halliday et al. (2001) along the 
major and minor axis of the galaxies revealed the presence of de- 
coupled kinematics on the radial scales of the NSDs. A counter- 
rotating central structure is present in NGC 4458, a cold rotating 
component in NGC 4478. 

(iii) The Lick/IDS absorption line-strength indices H/3, Mg h 
and (Fe) measured as a function of radius along the major and the 
minor axis of the galaxies allow us to estimate the age, metallicity 
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Figure 7. Top panels: The line strength of the H/3 and [MgFe]' indices at different positions in the galaxies NGC 4458 (left panel) and NGC 4478 (right 
panel). The grid shows the models of Thomas et al. (2003) for different metallicities [Z/B.] and ages. Bottom panels: The line strength of the Mg b and (Fe) 
indices at different positions in the galaxies NGC 4458 (left panel) and NGC 4478 (right panel). In all panels the circles show the major axis, the triangles 
the minor axis, the filled symbols the galaxy centres. The grid shows the models of Thomas et al. (2003) for different overabundances [a/Fe] and metallicities 
[Z/K], with age 12 Gyr. 



and overabundance of the stellar populations of the two galaxies 
using simple stellar population models. 

(iv) The radial variations of the derived quantities constrain the 
stellar population properties of the central discs. The NSD of NGC 
4458 has an estimated stellar mass of 3.4 x 10^ Mq and is similar to 
the decoupled cores of bright ellipticals, being as old as and richer 
in metals than the rest of the galaxy, with high and approximately 
constant overabundance. The cold central disc of NGC 4478 has 
a similar estimated stellar mass of 2.5 X 10'^ Mg, but IS younger, 
richer in metals and less overabundant than the main body of the 
galaxy. 

(v) The nearly solar q— element abundance of the central disc 
of NGC 4478 indicates a prolonged star formation history, typical 
of an tmdisturbed disc-like, gas-rich (possibly pre-enriched) struc- 
ture. 
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